%求解控制序列为a_max的相关参数
%t_f = x(1)
%编写一个计算 F 的函数，即 x 处方程的值
            %中间变量t1时刻的速度与位移
            F = @(x)[12+3*x(1)-12*(1 - 0.15 * sin(0.2 * x(1)));
                30+12*x(1)+0.5*3*x(1)*x(1)-(2*10-2*1.8*12*(1 - 0.15 * sin(0.2 * x(1)))+(100 + 12 * x(1) + 10 * cos(0.2 * x(1)) - 10))];
            
        %创建初始点 x0 
        x0 = [0];
        %设置选项以返回迭代显示。
        options = optimoptions('fsolve','Display','iter');